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Abstract 

Quantum adiabatic algorithm is a method of solving computational problems by evolving the 
ground state of a slowly varying Hamiltonian. The technique uses evolution of the ground state 
of a slowly varying Hamiltonian to reach the required output state. In some cases, such as the 
adiabatic versions of Grover's search algorithm and Deutsch-Jozsa algorithm, applying the global 
adiabatic evolution yields a complexity similar to their classical algorithms. However, using the 
local adiabatic evolution, the algorithms given by J. Roland and N. J. Cerf for Grover's search [ 
Phys. Rev. A. 65 042308(2002)] and by Saurya Das, Randy Kobes and Gabor Kunstatter for the 
Deutsch-Jozsa algorithm [Phys. Rev. A. 65, 062301 (2002)], yield a complexity of order vÏV (where 
N=2 n and n is the number of qubits). In this paper we report the experimental implementation 
of these local adiabatic evolution algorithms on a two qubit quantum information processor, by 
Nuclear Magnètic Resonance. 
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1. INTRODUCTION 



Quantum algorithms provide elegant opportunities to harness available 
quantum resources and perform certain computational tasks more efficiently 
than classical devices. The idea that a quantum computer could simulate 
the physical behavior of a quantum system as well as perform computation, 
attracted immediate attention l2| . The theory of such quantum computers 
is now well understood and several quantum algorithms like Deutsch-Jozsa 



n 



(DJ) algorithm |3], Grover's search algorithm 4], Shor's prime factoriza- 
tion algorithm Hogg's algorithm p ,Bernstein-Vazirani problem [7] and 
quantum counting |8| have been developed . All these algorithms start from 
a well-defmed initial state and perform computation by a sequence of re- 
versible lògic gates. After computation, the final state of the system gives 
the output. Various methods are being examined for building a quantum in- 
formation processing (QIP) device which is coherent and unitary J9|. Nuclear 
Magnètic Resonance has emerged as a leading candidate for implementation 



of various quantum computational prob 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 



ems on a small number of qubits 
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Quantum adiabatic algorithms provide an alternative method for comput- 
ing j3, í 



In this method the computation is done by evolving the system 
under a Hamiltonian for a given amount of time. Such algorithms start 
from a suitable input ground state and by evolution under a slowly time- 
varying Hamiltonian, reach the desired output state. Quantum adiabatic 



algorithms 



lems (27, 



lave been efficiently applied to solve various optimization prob- 



28 



29 



301 ]. Chuang et al. have demonstrated the implementation of 



a quantum adiabatic algorithm by solving the MAX-CUT problem on a 



three qubit system by NMR 



. In these algorithms, the condition for adi- 
abaticity is fulfilled globally by using only the minimum energy gap between 
the ground state and the first excited state for calculating the time of evolu- 
tion. This method of evolution is not efncient in some cases such as adiabatic 
Grover's search algorithm and adiabatic Deutsch-Jozsa algorithm as they re- 
sult in a complexity O(N) (N is the size of the data set), which is as good 
as their classical algorithms. However, these algorithms can be improved 
by application of local adiabatic evolution, where the adiabatic condition is 
fulfilled at each instant of time. This technique has been adopted theoret- 



ically by Roland and Cerf |37| for the adiabatic Grover's search algorithm 
and by S. Das et al. for adiabatic Deutsch-Jozsa algorithm |33| yielding 
a complexity 0(y/~N). Experimental implementation of adiabatic Grover's 
search algorithm based on the proposal of Roland and Cerf and adiabatic 
Deutsch-Jozsa algorithm of S. Das et al, is reported here. Section 2 contains 
an introduction to adiabatic algorithms. Section 3 discusses the adiabatic 
version of the Grover's search algorithm proposed by Roland and Cerf and 
its NMR implementation. Section 4 discusses the adiabatic Deutsch-Jozsa 
algorithm and its NMR implementation. Section 5 contains the experimental 
results, on a 2-qubit system, for both these algorithms. To the best of our 
knowledge this is the first experimental implementation of adiabatic Grover's 
search and adiabatic Deutsch-Jozsa algorithms. 



2. ADIABATIC ALGORITHM 



The adiabatic theorem of quantum mechanics states that when a system is 
evolved under a slowly time varying Hamiltonian, it stays in its instant aneous 



ground state 



28 



29 



34{ . This fact is used in solving certain còmput ational problems 
The problem to be solved is encoded in a final Hamiltonian 



(Hp), whose ground state is not easy to find. Adiabatic algorithms start with 
the ground state of a beginning Hamiltonian (Hb) which is easy to construct 
and whose ground state is also easy to prepare. The ground state of H B , 
which is a superposition of all the eigenstates of H^, is evolved under a time 
varying Hamiltonian H(s). H(s) is a linear interpolation of the beginning 
Hamiltonian Hb and the final Hamiltonian Hp such that 



H(s) = (l-s)H B + sH F , 



where < s < 1. 



(i) 



The parameter s = t/T tota i, where T tota i is the total time of evolution and t 
varies from to Ttotai- After evolution under the Hamiltonian H(s) for a time 
Tt tah the system is in the ground state of Hp with a probability (1 — £ 2 ) 2 , 
provided the evolution rate satisfies, 



max 

0<s<l 



dH(s) 
dt 



0;s 



9 



mm 



(2) 
The 



and the parameters of the algorithm are chosen to make e «Cl 
numerator in Eq. 2 is the transition amplitude between the ground state 
and the first excited state of H(s), and the denominator is the square of 
the smallest energy gap (g m i n ) between them. Ideally the time of evolution 
{Ttotai) must be infinite. However as long as the gap is finite, for any finite 
and positive e, the time of evolution can be finite. The time of evolution 



of the algorithm is determined by the minimum energy gap between the 
ground state and the first excited state. In the adiabatic case the time of 
evolution determines the complexity of the algorithms (that is how long it 
takes for the task to be completed), which can then be compared to the 
complexity of the discrete algorithms in classical and quantum paradigms. 
The time of evolution is measured in units of natural time scale associated 
with the system, T which is 0(h/E) where E is the fundamental energy scale 



associated with the physical system used to construct the states. 33] . 

In the actual implementation, the Hamiltonian H( s) is discretized into M+ 
1 steps as H(jj) where m goes from — > M J32|, |35|. Thus the time varying 
Hamiltonian H(s) goes from beginning Hamiltonian to final Hamiltonian in 
M+l steps. As the total number of steps increase, the evolution becomes 



more and more adiabatic 
given by 



The evolution operator for the m th step is 



U = e 



|(l-B)HB+BHr]Ai 



(3) 



where At = T/(M + 1). The total evolution is given by, 



M 



u = n u n 



(4) 



m=0 



Since, Hb and Hp do not commute in general, the evolution operator of Eq. 
3 is approximated to first order in At, by the use of the Trotter's formula 



32] as 



m 



(5) 



Thus in each step only a small evolution of the system from ground state of 
Hb towards the ground state of Hf takes place. 



3. GROVER'S SEARCH ALGORITHM 



Suppose we are given an unsorted database of N items and one of those 
items is marked. To search for the marked item classically, it would require 
on an average N/2 queries. However using quantum resources, the algorithm 

n i— 

prescribed by Grover |4j performs the same search with O(viV) queries. 
The algorithm starts with an equal superposition of states, representing the 
items, repeatedly flips the amplitude of the marked state (done by the oracle) 
followed by the flip of the amplitudes of all the states about the mean. The 
number of times this process is repeated determines the complexity of the 
algorithm and this scales with the size of the database as 0(\/N). 

In the adiabatic version, the system is evolved under a time dependent 
Hamiltonian which is a linear interpolation of H# and H^. As n qubits are 
used to label a database of size N (=2"), the resulting Hilbert space is of 
dimension N. The basis states in this space are \i) where i=0,- • • ,N. is 
chosen such that the ground state is a linear superposition of all the basis 
states. Therefore for a 2-qubit case, 



\^b) 



1(|00> + |01) + |10) + |11)) 



(6) 
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\iPb){iPb\, 
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(7) 



The Final Hamiltonian has the marked state \i/jf) as the ground state. 



h f = i- |<M(<M 



(8) 



The rate at which the interpolating Hamiltonian H(s) (given by Eq. [T]) 
changes from Hb to Hp depends on the condition, 

ds 



dt 



< £ 



9 2 {s) 



ds ' 



(9) 



Following Roland and Cerf j^Zj, t is obtained as a function of s as, 
1 N 



t 



arctan{\/ N — 1 (2s — 1)} + arctan\/ N — 1 



2e y/W^Ï 

Taking t' = et and on inverting the above function, s(t') is obtained as 



(10) 



4t') - \ 



{ ; -tan { — arctanV N — l) } + 1 

WW^Ï V N J 



(11) 



The plot of this function for N=4 (for a 2 qubit case) is given in Fig. 1. In 
the experiment the time of evolution is varied according to Eq. [TU It has 
been shown by Roland and Cerf |37| that with this adiabatic evolution, the 
complexity of the algorithm is O(VN). 



3.1. EXPERIMENTAL IMPLEMENTATION 



The NMR Hamiltonian for a weakly coupled two-spin system is 



(12) 



where lü\ and lü2 are Larmour freqüències and J\i the indirect spin-spin cou- 
pling. The beginning Hamiltonian for a 2-qubit Grover's algorithm as stated 
in Eq. UI written in terms of spin-half operators, is 



3 1 



(13) 



The identity term does not cause any evolution of the state and so it can be 
omitted, yielding the beginning Hamiltonian without the negative sign and 
the factor half as: 



7~Íb — Ixi + Ixi + 2/ x i7 X 2 



(14) 



The evolution under 7Ïb can be simulated by a free evolution under the 
Hamiltonian 7ï of Eq. [T2]between two 7r/2 pulses with appropriate phases. 



^(Iyt+Iyl) . p ÏHT m p -Íj(Iyl+Iy2) _ p Í{ L·ü l I xl+^2lx2 + ^J IxlIx2)T 



AH'T 



Let the state 1 00) be the marked state. The final Hamiltonian is, 

(\ 



(15) 



H 



|00) 
F 



I 








(16) 



\0 0/ 

In terms of spin operators the final Hamiltonian is, 



(17) 



The final Hamiltonian keeping the spin operator terms only and without the 
negative sign and the factor half is 



H 



|00) _ 



F 



hi + Iz2 + 21zllz2 



(18) 



Similarly the final Hamiltonian for other states being marked, in terms of 
the spin-half operators, is 

(19) 



1-èp 1 " 1 = I z i — I Z 2 — 2I Z \I Z 2 



H 



10) _ 



F 



111) 



W' = -Li 



Izl + Iz2 — 21zllz2 
Iz2 + 21zllz2 



(20) 
(21) 



The schematic representation of the experiment for the adiabatic Grover's 
algorithm in a two qubit system [consisting of a 1 H spin and a 13 C spin] is 
shown in Fig. 2a. The experiment is divided into three parts. The first 
part (preparation part) consists of preparation of pseudo-pure state (PPS) 
followed by equal superposition. The second part is the adiabatic evolu- 
tion, and the third part is the tomography of the resultant state. The pulse 
programme for the preparation of PPS and equal superposition is shown in 
Fig. 2b. The PPS is prepared by the method of spatial averaging |38j ] . After 
preparing PPS, equal superposition of states is obtained by application of the 
Hadamard gate on both the qubits. The Hadamard gate is implemented by 
(V/2 V - pulses, followed by tt x -pulses on both proton and carbon spins (Fig. 
2b) jl5| . The next stage consists of adiabatic evolution which has been car- 
ried out in the present work in 60 steps. Each step of the adiabatic evolution 
(Figs. 2c, 2d, 2e and 2f) consists of evolution under the final Hamiltonian 
for a time r sandwiched between two evolutions under the beginning Hamil- 
tonian for a time (T-r)/2. T is the total evolution time for one step and is 
equal to 1/ttJ. The value of r(= s x -M varies from to T takes place as 's' 
increases from to 1 according to Eq. HH m 60 steps. The pulse sequence 
for the beginning Hamiltonian is a free evolution of the system juxtaposed 
between two 7r/2 pulses with appropriate phases on each of the spins (the 
part marked as He in Figs. 2c-2f). The pulse sequence for the final Hamilto- 
nian depends on the marked state as stated in Eqs. ÍT^H2T1 If the state |00) is 
the marked state, then the pulse sequence for the implementation of the final 
Hamiltonian is a free evolution of the system under the NMR Hamiltonian 
juxtaposed between two tt pulses on each of the spins (Fig. 2c). Similarly, 



if the state 1 1 ) is marked the pulse sequence for the final Hamiltonian is a 
free evolution of the system between two tt pulses on the spin 1 (Fig. 2d), if 
the state 1 10) is marked then the pulse sequence is a free evolution between 
two 7r pulses on the spin 2 (Fig. 2e) and if the state 1 1 1) is marked, then 
the pulse sequence simulating the final Hamiltonian is just a free evolution 
of the system under the NMR Hamiltonian (Fig. 2f). 

The third stage of the experiment is the tomography of the final density 
matrix after the adiabatic evolution. The density matrix of a 2-spin system is 
a 4x4 matrix consisting of 6 independent off-diagonal complex elements (the 
remaining 6 are their complex conjugates), and the four diagonal elements 
which are the populations of the various levels. The diagonal elements are 
measured by 90° pulses on each qubit preceded by a gradient pulse. The 
six off-diagonal elements consist of four single quantum (SQ), one double 
quantum (DQ) and one zero quantum (ZQ) coherences. The real and the 
imaginary SQ, DQ and ZQ coherences in terms of the spin operators are; 



SQ r t eal 


— lix i 2>(lixljz) 5 


Q^imag 


= Iiy Í 2(liylj z ), 


BQ real 


'^'(lixljx liyljy) 


DQ ima9 


= ^(liyljx H - lixljy) 


ZQ real 




ZQ tma9 





where i^ j = 1,2 represents the qubits. Although the single quantum terms 
are directly observable, for proper scaling, all the off-diagonal elements are 



observed by a common protocol of two experiments; 



A : 



B : 




(23) 



(24) 



where 9 denotes the pulse angle, 0i, 02 the pulse phases and G z a gradient 
pulse. The first two pulses of the experiment A (depending on the pulse 
angle 9 and the pulse phases 4>\ an d 02 ) convert terms like l{ a + 2I ia Ijp 
into diagonal terms given by + 2I iz Ij z , where a and (3 denote the x, 
y, or z component of the spin operators of the first and the second qubit 
respectively. The gradient destroys all the transverse magnetization retaining 
only the longitudinal terms. The last pulse converts the retained longitudinal 
magnetization I iz + 2I iz Ij z into observable terms I ix + 2I ix Ij z . Thus the 
magnitude of Ii a +2I ia Ijp is mapped on to Ii x +2I ix Ij z which is then observed. 
In experiment B, a 7r-pulse is applied on the spin 'j' just before the tt/2 pulse 
on the spin 'i'. This creates the observable term I ix — 2I ix Ij z . The sum and 
difference of the two experiments yields 2I ia and 2I ia Ijp respectively. Six 
different experiments are needed to be performed to map the whole density 
matrix (real and imaginary). The various pulse angles and phases required 
during the experiment, and the resultant terms that are observed due to them 
are given in Table I. Experiments I and II yield the SQ, and experiments III- 
VI yield the ZQ and DQ coherences. 

4. DEUTSCH-JOZSA ALGORITHM 

The Deutsch-Jozsa Algorithm determines whether a binary function f(x), 



/(4r € {0, 1}») {0, 1} 



is Constant or Balanced [36]. A constant function implies that the function 
has the same value or 1 for all x. A balanced function implies that the 
function is for half the vàlues of x and 1 for the other half . For a two 
qubit case the constant and the balanced functions are given in Table II. 

In the adiabatic version of the Deutsch-Jozsa algorithm, the beginning 
Hamiltonian and its ground state, for a two qubit system, is given by Eq. 
|7] and Eq. |ü] respectively. The final Hamiltonian is given by Eq. El and the 
ground state of the final Hamiltonian for two qubits is of the form 



hM = a|00) + -^(|01) + |10> + |ll», 



(25) 



where 



1 



a= 4 
(3 2 = 1 - a 2 



l)/(oo) + (_i)/(oi) + (_i)/(io) + 



(26) 



From Eq. ESlit is seen that when a = 1 the function / is constant, and when 
a = then it is balanced. Thus a is chosen depending on whether the func- 
tion to be encoded in the final Hamiltonian is constant or balanced. Using 
Eqs. iïU7IISPo1 and [27)1 the matrix for the interpolating Hamiltonian(i7(s)) can 
be written as |33j |; 

í\ 1 1 l\ Í3a (A 
1111 
1111 



H(s) = I 



1 



s 
3 



0/3/3/3 
(3 (3 (3 
\0 (3 (3(3) 



(27) 



V i i i i/ 

S. Das et al. have shown that on evolution under the Hamiltonian H(s) 



takes the initial state to the solution state \i(f) |33|] . In the next section 



we describe an NMR implementation of the above algorithm. 



4.1. NMR IMPLEMENTATION 



The adiabatic Deutsch- Jozsa algorithm also, is implemented on the 2-qubit 
system. The beginning Hamiltonian in terms of the spin-half operators is the 
same as given in Eq. HU, and its implementation has been discussed in section 
3.1. 

The final Hamiltonian, obtained from Eqs. d El ESI and EBl for constant 
case (a=l) yields, 

(\ 



H C F = I 





\0 0/ 



(28) 



and for balanced case (a=0) yields, 



1 



H° F = I-- 

F 3 







(29) 



111 
111 

\o 1 1 1; 

The above final Hamiltonians in terms of spin-half operators can be written 
respectively as, 



3 1 



(30) 



and, 



3 , 1 

4" 3 



(31) 



As the identity does not cause any evolution of the state we consider only 
the spin operator terms. Thus the final Hamiltonian keeping only the spin 
operators (dropping the minus sign), for the constant case, can be written as 

(32) 



1 



Hp = -j\Jz\ + Iz2 + 2/21/22}, 



and for the balanced case as 



1 



T~(>F — ~q(IzI +/z2 + 2/21/22) + s(lxllx2 + Iylly2) 



_l_lr 1 1 r 



x2 



l(/a;l/z2 + Izllx2j 



(33) 



The signs of Eqs. HHE2]andES]are changed for consistency. Since the various 
terms in Eq. ESI do not commute, the evolution under this Hamiltonian 
would require a complex pulse sequence in NMR. However, we have found 
that by keeping only the diagonal terms in the Eq. ESI the pulse sequence 
simplifies considerably with the information regarding the balanced nature 
of the problem still encoded in it. This truncated final Hamiltonian for the 
balanced case is given by; 

(H b F ) trunc = -\{I z i + I z2 + 2I zl I z2 ) (34) 
o 

The opposite signs of Eq. and Eq. EU distinguish the constant and the 
balanced case. In the following we show that the balanced nature of the 
Deutsch-Jozsa problem is still encoded in (H b F ) trunc . Substituting a = and 
(3=1 and dropping the off-diagonal terms from the last part of Eq. EU , we 
obtain 

(\ 1 1 i\ ío 

1111 s 10 
1 1 1 1 3 1 
\1 1 1 1) \0 1) 



H\s) = I 



1 



(35) 



The eigenvalues of this Hamiltonian are: 





1 - 


Ao = 




6 - 




1 " 


Ai = 






6 - 


A 2 = 


A 3 



3 + 2s - x/9 + s(7s - 15) 



3 + 2s + x/9 + s{7s - 15) 



3 



(36) 

(37) 
(38) 



The vàlues of Ao, Ai, A2 and A3 as a function of 's' are plotted in Fig. 3. Ao is 
the ground state. As 's' increases from 0, Ao continues to be the ground state 
and becomes the ground state of the final Hamiltonian in the limit s — ► 1. 
The eigenvectors corresponding to Ao, Ai, A2 and A3 are respectively obtained 
as; 

1 " 

1 
1 

The final state to which the system converges after the evolution is 

/ \ 



( 3-s-2x/9-15a+7s 2 \ 
3(s-l) 

1 
1 
1 



V 



/ 3-s+2V9-15s+7s 2 \ 
3(*-l) 



I 



\ 



, v 2 - 



I 



-1 



V 1 / 



^3 



-1 
1 

V / 



lim vq = 

s— >1 



1 
1 

V 1 / 



(40) 



which is the desired output state. 

The energy gap between the ground state and the states corresponding 
to A2 and A3 goes to zero as s — > 1 as shown in Fig. 3. However, there 
is no transition from Ao to A2, A3 as the transition amplitude given by the 
numerator in Eq. [21 is zero in these cases. Therefore the transition amplitude 



from the ground state Ao to the next excited state Ai is relevant for calculation 
of s(t). The minimum energy gap between Ao and Ai, needed in Eq. El is 
obtained for s ~ 1 as seen in Fig. 3. Since the algorithm is implemented 
using local adiabatic evolutions we need to change s(t) such that the adiabatic 
condition 



ds <£ \g(s)\ 



dt -~\m 

as 



(41) 



is met at each time interval. Here g(s) is the energy gap between the ground 
state and the first excited state, given by |-\/9 — 15s + 7s 2 and \ {dH/ds) \ = 
Hp — Hb- The Hamiltonian is evolved at a rate that is a solution of 



ds ^ \g(s)\ 
dt £ \H F 



H 



B 



(42) 



On integrating Eq. 



we obtain í as a function of s. 
1 Us - 15 



t = 



+ k, 



(43) 



£2 v / 3v / 7s 2 - 15s + 9 
where the constant of integration k = to obey s = at t = 0. Inverting 
this function we obtain s(t') as 



s(t') 



3 

Ü 



^225 + 24í (55^ - 183t + 60V3t 2 - 18í 3 ) 



3 + 20V3t - I2t 2 



(44) 



where t' is et. The plot of s as a function of t! is shown in Fig. 4. From Figs. 
3 and 4 it is seen that the rate of change of s (and hence of the Hamiltonian) 
is fast when the energy gap between Ao and Ai is large, and slow when the 
gap is small. In practice the time of evolution for Hb and Hp is given by 
(1 — s) X T and s x T respectively, where T is 1/ttJ and s is varied from 
to 1 according to Eq. H31 In our implementation, the t interval for which s 



varies from to 1 is divided in 80 equal steps, and the corresponding vàlues 
of s for each step (calculated from Eq. 48) are substituted in the evolution 
time of Hb and Hp. 

On integrating Eq. 46 from s=0 to s=l, we get the total time of evolution 

= (45) 

T^totai is given in the units of T which is the time scale associated with the 



physical system used |33| . The time scale associated with evolution under the 
NMR Hamiltonian is ~ 10~ 3 s. The total time of evolution of the experiment 
(Ttotal) is given by 80 xT, where T is the time for one step (see Fig. 5b). For 
the choice e ~ 10 -2 , T ~ 60 x 10~ 3 s in our case. 

4.2. EXPERIMENTAL IMPLEMENTATION 

The experimental implementation of adiabatic Deutsch-Jozsa algorithm 
on a 2-qubit system [consisting of a 1 H spin and a 13 C spin] also consists 
of three parts namely preparation, adiabatic evolution and tomography of 
the final density matrix. The preparation of the pseudo pure state (PPS) 
and making of equal superposing of states as well as the tomography of the 
final states has already been discussed in section 3. So we only describe the 
method of implementation of the final Hamiltonian for the Deutsch-Jozsa 
algorithm. 

The pulse sequence for the implementation of the constant case final 
Hamiltonian (Tip) is given in Fig. 5b. The beginning Hamiltonian is im- 
plemented by a free evolution juxtaposed between tï/2 pulses with required 
phases (Fig. 5b). The implementation of the final Hamiltonian for the con- 



stant case is a free evolution under the NMR Hamiltonian of Eq. 14, juxta- 
posed between two 7r-pulses as shown in Fig. 5b. In the balanced case the 
implementation of the beginning Hamiltonian is same as in Fig. 5b. How- 
ever, the implementation of the final Hamiltonian (7i b F ) trunc is done in two 
parts [Fig. 5c]. The first part is a free evolution under the Hamiltonian 
given in Eq. 14 [Tf period in Fig. 5c]. The operator corresponding to such 
an evolution for time r will be of the form; 

e ÍTrJ(-I zl -I Z 2+2I zl I z2 )T {^Q) 

In the second evolution of 2r, the chemical shifts are refocused so that the 
system evolves only under its scalar coupling Hamiltonian 2ttJI z \I Z 2. Just 
before and after the evolution 7r-pulses with appropriate phases are put on 
each of the spins to flip the sign of the corresponding spin operator [Tj period 
in Fig. 5c]. The operator for the sequence of two pulses with an intermediate 
evolution for 2r is of the form 

e -i(Ixi)n . e ÍTrJ(2I zl I z2 )2T _ gí(Ixi)7T _ ç-mJ{2I zl I z2 )2r {^7) 

As these two evolutions given in Eq. EE1 and Eq. EU commute, the effective 
evolution for the 3r period is: 

e ÍTrJ(-I zl -I z2 +2I z iI z2 )T _ ç-íitJ(2I z1 I z2 )2t _ ^J{-I z \-1 z2 -21 z1 I z2 )t {^8) 

Thus the evolution during Tj cancels the J-evolution during Tf and adds 
a minus sign to it, yielding the effective Hamiltonian of Eq. HH1 and an 
effective evolution time of r. An evolution time of r=l/7rJ implements the 
full Hamiltonian of Eq. 39 as required for adiabatic evolution. Overall the 
cycle time for each step for the balanced case is increased to T+2r. 



5. EXPERIMENTAL RESULTS 



The experiments have been carried out using carbon-13 labeled chloroform 
( 13 CHCi3) where the two spins l· l·l and 13 C form the two qubit system. The 
proton spin represents the first qubit and carbon-13 the second. The sample 
of 13 CHCl3 was dissolved in the solvent CDCI3 and the experiments were 
performed at room temperature in a magnètic field of 11.2 Tesla. At this field 
the X H resonance frequency is 500.13 MHz and the 13 C resonance frequency is 
125.76 MHz. During the entire experiment, the transmitter freqüències of l· l·l 
and 13 C are set at a value J/2 away from resonance to achieve the condition 
W\ = u<2 = 7rJ. The equilibrium spectra of the two qubits are shown in Fig. 
6a, and the spectrum corresponding to |00) PPS is shown in Fig. 6b. To 
quantify the experimental result we calculate the average absolute deviation 
42] of each element of the experimentally obtained density matrix from each 
element of the theoretically predicted density matrix given by, 



where N=2 n (n being the number of qubits), xfj is element of the 

theoretically predicted density matrix and xfj is (i,j) th element of the ex- 
perimentally obtained density matrix. 

5.1. GROVER'S SEARCH ALGORITHM 




(49) 



The experimental spectra corresponding to the implementation of Grover's 
search algorithm on the above two qubit system are given in Fig. 7. the 
spectra given in Figs. 7a(i-iv) contain the reading of populations after re- 



spectively searching states |00), |01),|10) and |11). The population spectra 
are obtained by application of a gradient followed by a tt/2 pulse. Depending 
on the final state, the population spectra consist of one single spectral line 
for each spin. These correspond to, |00) — ► 1 01) and |00) — ► |10}transition 
when the searched state is |00) (Fig. 7a-i); |01) -> |00) and |01) |11) 
when the search state 1 01 ) (Fig. 7a-ii); 1 10) — > |00) and 1 10) — > |11) when the 
search state 1 10) (Fig. 7a-iii); |11) — > |01) and 1 11) — > 1 10) when the search 
state 1 1 1 ) (7a- iv). The coherence spectra in Fig. 7b have been obtained 
by observing the searched state without application of any r.f. pulses. The 
absence of any signal in the spectra confirms that there is no single quantum 
coherences after the search. To check for the absence of zero quantum and 
double quantum coherences as well, the entire density matrix has been tomo- 
graphed. Fig 8a shows the theoretical and the experimental density matrices 
after the adiabatic evolution, when state 1 00) has been searched. The mean 
deviation of the experimentally obtained density matrix from the theoreti- 
cally predicted one (calculated using Eq. HHJ) is 2.49%. Similarly Figs. 8b, 
8c and 8d contain the theoretically predicted and experimentally obtained 
density matrices when the states |01), 1 10) and 1 1 1) have been searched. The 
mean deviation of the experimental density matrices from their theoretically 
predicted counterparts are 1.92%, 1.89% and 1.97% respectively. 

5.2. DEUTSCH-JOZSA ALGORITHM 

5.2.1 Constant case 

For the constant case (Eq. EBD, the state expected after the evolution (using 
the pulse sequence given in Fig. 5b) is |00). The density matrix consists 



of population in 1 00) state and no coherences. The spectrum correspond- 
ing to the population for such a state, obtained by application of a gradient 
followed by tt/2 pulses on each of the spins, consists of one single quantum 
coherence in each spin ('Population spectrum' in Fig. 9a). The spectrum for 
coherence, observed without application of any pulses on any of the spins, 
has a near absence of any signal ('Coherence spectrum' in Fig. 9a). Further 
confirmation of the final state is done by the tomography of the complete 
density matrix. The Fig. 10 shows the tomography of the experimental and 
theoretically predicted density matrices of the final state for the constant 
case. The mean deviation of the experimental density matrix from the the- 
oretical one is 5.28% 
5.2.2 Balanced case 

For the balanced case [Eq. 31, a=0 and (3=1], the state expected after the 
evolution (using the pulse sequence of Fig. 5c) is -j= ( 1 01) + 1 10) + |11)). The 
theoretical density matrix of the final state is given in Fig 11 (a). This state 
theoretically has three diagonal elements, one SQ coherence of each qubit 
and a ZQ coherence between the two qubits, all of equal intensity. This 
state is confirmed by the spectra shown in Fig. 9b and the density matrix 
in Fig. ll(b). The mean deviation of the experimentally obtained density 
matrix from the theoretically predicted one is 17.2%. It is seen that in the 
density matrix obtained from experiment, the SQ coherence of 13 C (second 
qubit) and the ZQ coherence between 13 C and 1 H have significantly reduced 
intensity, compared to the theoretically expected vàlues. 

There are three sources of error in adiabatic algorithms. e gives a mea- 
sure of the first source of error. Theoretically the total time of evolution in 



adiabatic algorithms should be infinite. However, in practice the evolution is 
terminated once the state is supposed to have been reached with sufficiently 
high probability given by (1 — e 2 ) 2 which in our case (for e = 10~ 2 )is obtained 
to be 99.98%. The second source of error is due to neglect of 0(Aí 3 ) terms 
in the Trotter's Formula (Eq. EJ). The maximum error introduced due to this 
is « 0.92 % which can be safely neglected. 

The third source of error is due to decoherence effects arising from the 
interaction of the spins with their surroundings. To study decoherence, the 
relaxation times Ti and T2 of 1 H and 13 C were measured. The T2 for SQ 
coherences were measured by CPMG sequence. For the measurement of ZQ 
and DQ coherence decay rate, the term li x ^2x was created and its relaxation 
rate was measured by CPMG sequence. The T2 of SQ coherence of l· l·l was 
found to be 3.4 s and for 13 C it was found to be 0.29 s. The decay rate of I\ x ï-2x 
term was found to be 0.19 s. The Ti for l· l·l and 13 C measured from the initial 
part of the inversion recovery experiment was found to be 21 s for 1 H and 
16s for 13 C. Using these measured vàlues of Ti and T2 the simulation for the 
balanced case was repeated including relaxation using Bloch's equations 
Significant decay of the carbon coherences was observed. The mean deviation 
of the of the experimental density matrix form the theoretical density matrix 
including relaxation is found to be 8.0%. 

The observed mean deviation between the theoretically expected and the 
experimentally obtained density matrices for the Grover's search and the 
constant case of the Deutsch-Jozsa are small (< 2% and < 6% respectively) 
while that for the balanced case of the Deutsch-Jozsa is large (~ 17%). In the 
first two cases, the results are encoded in the diagonal elements of the density 



matrix, which are attenuated by the spin lattice relaxation, the times for 
which are large (> 16 sec). On the other hand, in the balanced Deutsch-Jozsa 
case, there are off-diagonal elements as well which are attenuated by spin- 
spin relaxation, the times for which are small (< 4 sec for l· l·l and < 0.3 sec for 
13 C). The decoherence times thus have a large effect in this case. A correction 
for the decoherence has improved the mean deviation considerably (reduced 
to ~ 8%), confirming the succesful implementation of these algorithms. 

6. CONCLUSION 

In this paper we have demonstrated the experimental implementation of 
Grover's search and Deutsch-Jozsa algorithms by using local adiabatic evolu- 
tion in a two-qubit quantum computer by nuclear magnètic resonance tech- 
nique. We have suggested a different Hamiltonian for the adiabatic Deutsch- 
Jozsa algorithm which is diagonal in the computational basis and hence easier 
to implement by NMR. To the best of our knowledge this is the first exper- 
imental implementation of these two algorithms by adiabatic evolution.We 
believe that this work will provide Ímpetus to solving other problems by adi- 
abatic evolution. 
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TABLE CAPTIONS 



Table I: The phases 4>\ an d 02 and the flip angle of the 6 pulses of the exper- 
iments A and B (see Eq. 26 and 27). a and (3 denotes the various 
components of the spin operator terms of the first and the second qubit 
whose magnitudes were determined by that particular experiment. 

Table II: The two constant and the six balanced functions for the 2-bit Deutsch- 
Jozsa algorithm. 
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FIGURE CAPTIONS 



Fig 1: The plot of s as a function of t' for the 2-qubit adiabatic Grover's search 
algorithm (ref. Eq. 13). s(t') is plotted for that interval of t' in which 
's' goes from — > 1. 

Fig 2: Pulse sequence for the implementation of adiabatic Grover's search al- 
gorithm. The narrow filled pulses, which donot have any angle specified 
on them, represents 90° pulses while the broad unfilled pulses represents 
180° pulses. The phases (X and Y represents -X and -Y phases respec- 
tively) are specified on each pulse. (a) Schematic representation of the 
pulse programme. The part 'preparation' consists of creation of PPS, 
followed by equal superposition. The second part is adiabatic evolu- 
tion done in 60 iterations and the third part is tomography of the final 
density matrix. (b) Pulse sequence for preparation of PPS and equal 
superposition. (c) Pulse sequence for the adiabatic evolution when the 
|00) is been searched. In this sequence the system is evolved under H# 
for time T-r and under for a time r. This is repeated 60 times for 
various r varying it from — > T as s is varied from — > 1 according to 
the Eq. 13 in equal intervals of t'. The total time in each iteration is T 
(= 1/VJ) which is 1.52 ms. The pulse sequences (d),(e) and (f) when 
Hp is encoded to search for the states |01), 1 10) and |11) respectively. 

Fig 3: The eigenvalues of H b (s) for the Deutsch-Jozsa algorithm (Eqs. 41-43) 
plotted as a function of parameter s. Xq is the ground state. Ài,À2 and 
A3 are the excited states. A2 and A3 are degenerate for all vàlues of s. 
Xq approaches A2 and A3 as s— >1. Ai changes marginally as a function 



of s. 

4: Plot of the parameter s as a function of t' for the Deutsch-Jozsa algo- 
rithm (Eq. 48). s is for t'=0 and s is 1 for t' = -^.This shows that 
s changes rapidly at the beginning when |Ao — Ai| is large (Fig. 1) and 
later it changes slowly as |Ao — Ai| becomes small. 

5: Pulse programme for the implementation of the adiabatic Deutsch-Jozsa 
algorithm. The narrow filled pulses which do not have any angle spec- 
ified on them are 90° pulses and all the broad unfilled pulses are 180° 
pulses. The phase of the pulses are specified on each of them. The fre- 
quency offset is set at a value J/2 during all the evolutions. (a) Block 
diagram representation of the pulse programme. The preparation se- 
quence has already been explained in Fig. 2. The adiabatic evolution 
is shown in (b) and (c). The measurement process is the tomography 
of the final density matrix which is explained in the text. (b) Pulse se- 
quence of the adiabatic evolution under the interpolating Hamiltonian 
H(s) for the constant case. Hb represents the beginning Hamiltonian 
and the pulse sequence implements the evolution as given in Eq. HH1 
Hp represents the final Hamiltonian when the constant case in encoded 
in it. T is the effective time of evolution in each cycle and r goes from 
— > T slowly in 80 steps. (c) The pulse sequence of H(s) for the bal- 
anced case. Hp 1 represents the pulse sequence for the implementation 
of the final Hamiltonian when the balanced case is encoded in it. Dur- 
ing the period Tf, the evolution takes place under the free Hamiltonian 
given by Eq. [T2J During the period Tj, the evolution takes place under 



the J-coupling Hamiltonian given in Eq. USI The 7r-pulse between Tf 
and Tj and at the end of the Tj restores the correct sign of the coupling 
Hamiltonian such that the total evolution for 3r is given by and the 
effective time of evolution for Hp 1 is just r. The time r is incremented 
in 80 steps from — > T, where T (=l/7rJ) is 1.5 ms in our experiments. 

6: (a) Equilibrium spectra of 13 CHCl3. The small line in the middle of 
the 1 H spectrum is due to proton of unlabeled chloroform and the three 
small equal intensity lines in the 13 C spectrum are due to the J-coupling 
of the deuteron with the natural abundant 13 C in the solvent CDCI3 
(b)The spectra corresponding to |00) PPS. A single line with positive 
intensity on each of the spins confirm that only |00) level is populated. 

7: Results of Grover's search algorithm. (a) The population spectra of 
the final state after the search have been performed. The populations 
have been observed by a applying a gradient followed by a tt/2 pulse. 
Figs. 7a(i - iv) contain the population spectra respectively when the 
states |00), |01), 1 10) and 1 11) have been searched. (b) The spectra of 
coherences (obtained by observing without the application of any pulse) 
of the final states after the adiabatic evolution. Figs. 7b (i - iv) contain 
the coherence spectra when the search states are |00), |01), 1 10) and 1 1 1 ) 
respectively. 

8: The tomography of the real and imaginary parts of theoretically ex- 
pected and experimentally obtained density matrices for the search 
states in Grover's algorithm. (a) |01) (b) |01) (c) |01) and (d) |01). 
The density matrices consist of just a real term on the diagonal corre- 



sponding to the population of the state that has been searched. 

; 9: (a) Spectra of the population and coherence of 1 H and 13 C after the 
implementation of the constant case of the Deutsch-Jozsa algorithm. 
As the final state for the constant case is |00), the 'Population' spec- 
trum consists of one single quantum coherence for each spin and the 
'Coherence' spectrum contains no signal [see text], (b) The spectra of 
population and coherence for the balanced case of the Deutsch-Jozsa 
algorithm. The expected state is -^(|01) + 1 10) + |11)). As the popula- 
tion of 1 01} , 1 10) and 1 1 1) states are the same, the spectra consists one 
SQ coherence for each of the spins. The coherence spectrum consists of 
another observable SQ coherence for each spin [see text] . The observed 
intensity of the SQ coherences in the 'Population spectrum' is nearly 
half compared to those in the 'Coherence spectrum' according to the 
expectations. 

10: The tomography of the real and the imaginary parts of (a) theoreti- 
cally expected and (b) experimentally obtained density matrices after 
the implementation of the constant case of the adiabatic Deutsch-Jozsa 
algorithm. As the theoretically expected state is |00), the density ma- 
trices contains just one diagonal term which is the population of the 
|00) state. 

11: The tomography of the real and the imaginary parts of the (a) theoreti- 
cally predicted and (b) experimentally obtained density matrices of the 
final state for the balanced case of adiabatic Deutsch-Jozsa algorithm. 
1 H and 13 C are taken as the first and second qubit respectively. The 



theoretically predicted final state is ^=(|01)+|10)+|11)). Thereforethe 
density matrix contains three diagonal elements corresponding to the 
populations of |01), 1 10) and |11) states, SQ coherences corresponding 
to l· l·l and 13 C, and one ZQ coherence between the two qubits, all of 
equal intensity. (c) The real and the imaginary parts of the theoreti- 
cally calculated density matrix of the final state for the balanced case of 
the Deutsch-Jozsa algorithm after the inclusion of the relaxation effects. 
The decay rates used for the single quantum coherences are 3.4s for l· l·l 
and 0.29 s for 13 C. The decay rate of the zero quantum and double 
quantum coherences of this hetronuclear two spin system used is 0.19 s. 
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